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^ . Abstract 
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We define the chemical potential as the Lagrange multiplier of 

■ the baryon charge operator in the transfer matrix formalism of QCD 

^ . on a lattice. Transforming the partition function into an euclidean 

QJ \ path integral we get the Hasenfratz-Karsh action both for Wilson and 

' Kogut-Susskind fermions. In the latter case the chemical potential in 
the spin-diagonal basis is half that in the flavour basis. Some open 

* I 

^ , problems in the spin-diagonal basis are pointed out. 
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1 Introduction 



Hasenfratz and Karsh [|l| noticed that if in the euchdean path integral of 
relativistic field theories regularized on a lattice the chemical potential is 
coupled to the Fermi fields as in classical systems, specific counterterms are 
needed which would be extremely unconvenient in numerical simulations. 
Guided by the analogy between chemical potential and gauge fields, they 
exponentiated the chemical potential in the same way as the gauge fields 
are exponentiated in the link variables of Wilson regularization. They also 
showed that this definition does not require any counterterms (in the free 
theory), which makes it a good candidate for numerical simulations. 

From the point of view of renormalization the way the chemical potential 
is introduced is obviously not unique so that one is free to look for the 
form which is most convenient for a specific purpose. But Gavai gave an 
argument whereby in any case the chemical potential cannot couple linearly 
to the fermion fields in the path integral. 

There are a few points which in our view should be reexamined for a fully 
satisfactory settlement of the issue, the more so in view of the long standing 
difficulties met in numerical simulations 0]. 

For instance, the heuristic argument does not apply with a noncompact 
regularization ||^, where the gauge fields still couple linearly to the Fermi 
fields. Should one in such choose a different coupling of the chemical 

potential? Moreover in statistical mechanics the chemical potential is the 
Lagrange multiplier of a conserved operator, a feature which is lost if it is 
exponentiated. 

A last point to be clarified concerns the existence of the transfer matrix 
for the different fermion regularizations. This is of course a prerequisite for 
the lattice formulation to be correct, but its role becomes apparent in the 
present approach. We start from the observation that in classical statistical 
mechanics the partition function is the trace of the exponential of (minus) 
the hamiltonian. Now in the corresponding formulation of QCD the partition 
function is the trace of the transfer matrix, whose logarithm defines the 
hamiltonian. It is in this formulation that the chemical potential should 
behave as a Lagrange multiplier. We adopt this definition and show how the 
manipulations necessary to transform the trace into a path integral change 
the coupling of the chemical potential to the Fermi fields from linear to 
exponential. Following this prescription we find in fact exactly the result of 
Hasenfratz and Karsh both for Wilson and Kogut-Susskind fermions, with 
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the qualification discussed below. In the latter case the value of the chemical 
potential in the spin-diagonal basis is half that in the flavour basis, namely 
their ratio is the inverse of that of the cutoffs. 

In our derivation we assume a given form for the transfer matrix. Now for 
Wilson fermions the transfer matrix was constructed by Liischer in terms 
of quark- ant iquark creation-annihilation operators, whereby the hamiltonian 
and the baryonic charge find the proper definition. The proof of reflection 
positivity was given in the gauge Uq = 1. This gauge fixing is too strong, 
but Menotti and Pehssetto have extended the proof to the gauge [/q ~ 1, 
where [/q = 1 with the exception of a single time slice. 

For the case of Kogut-Susskind we have not found in the literature a 
transfer matrix with the same desirable features. In the only paper on the 
subject we are aware of, the transfer matrix is constructed in the spin- 
diagonal basis, where it is found necessary to introduce auxiliary fields whose 
relation to quark-antiquarks is not transparent. Another feature, possibly 
related to the above, might be a difficulty in the way of the application we 
are interested in: the transfer matrix is linear in the creation-annihilation 
operators, rather than exponential, so that the hamiltonian contains the 
logarithm of these operators. 

We have therefore tackled the problem in the flavour basis, and we have 
found that the formalism developed by Liischer ^ for Wilson fermions can 
easily be extended to the case of Kogut-Susskind, provided the elementary 
gauge variables are attached to the links of the blocks rather than of the 
lattice. This excludes the possibility of getting the action in the spin-diagonal 
basis from that in the flavour basis by a simple change of variables. If instead 
the elementary gauge variables are associated to the links of the lattice, we 
can transform to the spin-diagonal basis, but we meet with a formal difficulty 
and a practical complication. In fact the proof by Menotti and Pelissetto 
cannot be adapted straightforwardly, because their gauge fixing has in such 
a case a nontrivial Faddeev-Popov determinant (as explained in detail below) 
which must be taken into account to establish the reflection positivity of the 
transfer matrix. This is the formal difficulty. The practical complication is 
the following. We have constructed the transfer matrix in the flavour basis 
with the "minimal" gauge coupling, the one where the link variables join 
the quark fields along the shortest path. Now it is well known that a gauge 



coupling simple in one basis becomes very complicated in the other and 
in fact we find an awkward result in the spin-diagonal basis. 

Both shortcomings can hopefully be overcomed, but the construction of 
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the transfer matrix with Kogut-Susskind fermions in the spin-diagonal basis 
requires further study. 

Two final observations. The first is that no claim of uniqueness is made 
about the coupling of the chemical potential in the path integral formulation 
for a given transfer matrix. The question of uniqueness will be discussed in 
separate paper. 

The second is that the procedure just outlined appears also appropriate 
in the study of a system in the canonical ensemble, where since one has to 
integrate over the (imaginary) chemical potential [|12|, the exact way this is 
introduced can be of practical importance. 

The paper is organized in the following way. We first show how the present 
prescription reproduces the Hasenfratz-Karsh result for Wilson fermions, by 
using Liischer's construction of the transfer matrix. We then apply the 
same procedure to Kogut-Susskind fermions, adopting the notations of Mont- 
vay and Miinster |jlO|], and using the transfer matrix especially devised for 



the present application 0. 

In this paper the pure gauge part of the transfer matrix and the parti- 
tion function will be omitted because it does not play any explicit role (but 
with the Kogut-Susskind fermions in the spin-diagonal basis) and the lattice 
spacing is set equal to 1. 



2 Transfer matrix and chemical potential 

We start from the definiton of the grand canonical partition function at finite 
baryon density according to the ordered product 

Z = Tr|exp(|g^)n^.M}, (1) 

where the temperature T is the inverse of the number A^^o of links in the 
temporal direction, is the baryon charge operator and Tq is the quark 
transfer matrix. The chemical potential /i is determined by the condition 
that the expectation value of the baryon number be 

qs = Z-^TrlQBeMy^NoQB)X{%{n^)\ . (2) 

I no J 
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The transfer matrix Tq can be written 0, ^ in terms of an auxihary operator 
Tq^rio) 

%{no) = J-'Tj(no) T,K + s), (3) 

where JT" is a function of the gauge fields which will be defined later and 
s = ±1 for Wilson/Kogut-Susskind fermions respectively. This difference in 
sign does not reflect any intrinsic difference, but is only due to the different 
conventions adopted by Liischer and Montvay-Miinster, which we maintain 
for easy reference. As a consequence of these conventions the quarks prop- 
agate backwards in time with the Wilson definition and forwards with the 
Kogut-Susskind one. 

Tq is defined in terms of quark-antiquark creation-annihilation operators 
x^, y^, X, y acting in a Fock space. It depends on the time coordinate uq only 
throu the dependence on it of the gauge fields. In fact the creation and 
annihilation operators do not depend on no. They depend on the spatial 
coordinates n of the sites or, in the case of Kogut-Susskind in the flavour 
basis, of the blocks, and on Dirac, flavour and color indices, a, /, c sometimes 
comprehensively represented by /. 

In the transfer matrix formalism most often one has to do with quanti- 
ties at a given (euclidean) time tiq. For this reason we adopt a summation 
convention over spatial coordinates and intrinsic indices at fixed time. So for 
instance we will write 

g(mo)Q(mo,no)g(no) = ^m,/('^o)Qm,/;n,j("^0, '^o)'?n,j("'o), (4) 

m,n,/, J 

where Q is the quark matrix and q the quark field. Q is a function of time 
dependent link operators U^{nQ) which have the standard Wilson variables 
U^j,{mQ,m) as spatial matrix elements 

(f/^(mo))^_„ = 5^,JJi,{rnQ, m). (5) 

In this notation the quark action Sq, the baryonic charge Qb and the auxiliary 
operator Tq{no) can be written 

Sq= Y Qi^o)Qi'mo,no)qino) (6) 

mo, no 

Qb = x^x - (7) 
^j(^o) = exp — x^M(no) X — y^M(no) y ex.Y> [y N {tiq) x] . (8) 
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Their form is the same for Wilson and Kogut-Susskind fermions but the 
matrices Q, M and N are different in the two cases and will be specified 
later. The expression of Tg(no) is valid in the gauge ?7o = 1 which is not 
admissible. In the construction of the path integral formulation of QCD at 
finite baryon density we do not need to fix the gauge, but to lighten the 
formalism we will nevertheless put f/g = 1 and we will reinstate Uq in the 
final result. The reader can check that keeping Uq in the intermediate steps 
one arrives at the same result, provided some care is exercised: for instance 
when Uq ^ \ the expression x"^M(no) x + y^MiriQ) y appearing in Eq. (H) 
changes and does not commute with Qb any longer. We anticipate that N 
is hermitean and also M is hermitean in the gauge Uq = 1. 

Now we transform the trace into a Berezin integral. To this end we 
introduce between the factors in Eq. (|I]) the identity 

1 = J[dx~^dxdy^dy]ex.p{—x^x — y^y)\xy><xy\, (9) 

where the basis vectors 

\xy >= I exp(— xx^ — yy^) > (10) 

are coherent states and the x^,x,y'^,y are Grassmann variables. They are 
labeled by the time slice where the unit operator is introduced. For the 
other indices they are subject to the same convention as the creation and 
annihilation operators . 

We then use the following equations p. First 

< x\ exp (^x^Mctj \x' >= exp (^x^e^^ x'^ . (11) 

In particular 

< X, y\ exp {^iQb) y' >= exp (^e'^x^x' + e^^y^y'^ , (12) 

which shows how the chemical potential gets exponentiated going from the 
trace of the transfer matrix to the euclidean path integral. 
Second, for arbitrary matrices M and 

< x\ exp(x^Ma;) exp(i;^A^a;)|x' >= exp (^x^ e^^ e'^ x'^ . (13) 

Finally, if i? = i?(a;^) and C = C{x) are operators which depend on x^,x 
only 

< x\BMC\x' >= 5(x+) < x\M\x' > C{x'). (14) 
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Using these equations we can evaluate the kernel of the transfer matrix 

< x{no),y{no)\fq{no)\x{no + s),y{no + s) >= J^'^ exp (a;+(no)iV(no)|/+(no)) 

• exp (^x^{no) exp (— M(no)) exp (— M(no + s)) exp(/i) x{no + s)^ 

• exp (y^{no) exp (-M(no)) exp (-M(no + s)) exp(-;u) y(no + s)) 

• exp (y(no + s)N{no + s)x{no + s)) . (15) 

Collecting all the pieces we arrive at the euchdean path integral form of 
the partition function 

Z = J-^ J [dx^dx dy^dy] exp S' (16) 

with the action 

S' = exp{Noij)x+{l)x{l + s) + cxp{-Noi^)y+{l)y{l + s) 

+ ^ |(1 - Sno,i)x^{no) exp(-M(no)) exp(-M(no + s)) x{no + s) 

no 

+(1 - Sno,i)y^ino) exp(-M(no)) exp(-M(no + s)) y(no + s) 
+x^{no)N{no) y^{no) + y{no + s)N{no + s) x{no + s) 
-x+(no)x(no) - y+(no)y(no)} . (17) 

The chemical potential appears only on the first site. We can arrive at a 
uniform distribution by the change of variables 

x{no) exp {sa (no)) X (no), x^{no) cxp{— sa{no))x^ (no) , 

y{no) cxp{-sa{no))y{no), y^ (no) ^ exp{sa{no))y^ (no) (18) 

where ^ 

a{no)^a{l) + {no-l)fx--{l + s)Noi^, no > 1, (19) 

with a{l) arbitrary and a{no) = a{no + Nq) to respect the antiperiodicity in 
time of the quark field. We thus get 

yS' = — > ^ exp{iJ,)x'^{no) exp(— M(no)) exp(— M(no + s)) x{no + s) 

no 

+ exp(-/i)|/+(no)exp(-M(no)) exp(-M(no + s))y{no + s) 
+x^{no)N{no) y^{no) + y{no + s)N{no + s) x{no + s) 
-x+(no)a;(no) - y+(no)|/(no). (20) 
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3 Wilson fermions 

We set s = 1 in Eq. ( pOD and assume 

M(no) = -\n(^{2K)^ B'^rio] 
N{no) = 2KB{noyh{no)B{no)-K (21) 
where K is the hopping parameter and 

3 

B{no) = 1- Kj2(u,ino)T^^^ + Tj-^U+ino: 

3 

2 



cK) = lj:^crAu,{no)T^^^ -T^-^U^{no)). (22) 



We have introduced the translation operators T^^"* with matrix elements 

where the unit vectors have matrix elements 

(e,). = 5,,.. (24) 
Next we define the Dirac spinors q by the transformation 

X = B^P^'^q, y = B^Pt^q 

q = gt^o (25) 

where 

^0^^^ = ^ (1 ± 7o) . (26) 

The jacobian of this transformation is the function J' introduced in Eq. (^. 
The charge conjugation transformation is 

q = C-^q' 

q = -C^q' (27) 
where, with Liisher's conventions 

C = 7o72. (28) 
8 



It exchanges x and y according to 

X = 72?/' 

y = 12 x', (29) 

justifying our definition (|^) of the baryon charge. Now we reinstate Uq. The 
partition function takes the form 

Z = J [dqdq] exp 5*,, (30) 

with the action 

S, = ^{exp(/x)2irg(no)Po^+^f/o(^o)g(^o + l) 

no 

+ exp(-/i)2irg(no + l)Pt^Ut\no)q{no) 

+g(no) [2KCino) - 5(no)] g(no)} (31) 

where 

C{no) = ^ E 7. {UAno)Tj^^ - ^^"^f^/K)) • (32) 

This is the Hasenfratz-Karsch action for Wilson fermions with Wilson pa- 
rameter r = 1. Notice the "plus" sign in the exponential of the action, to 
comply with Liischer's convention. 



4 Kogut-Susskind fermions in the flavour ba- 
sis 

The Kogut-Susskind fermions can be defined in the flavour as well as in the 
spin-diagonal basis. If we want to be able to transform from one basis to 
the other, the gauge fields must be defined on the links of the lattice. But 
if we want to stay in the flavour basis, as we will do in this Section, we can 
associate the gauge fields to the links of the blocks and forget the lattice. 
We set s = —1 in Eq. (^) and assume ^ 



M(no) = 



3 

(+) 



-^i^^^i"^f^/K))l + ^l®l + 75® t,to, (33) 
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where m is the quark mass parameter and K the hopping parameter. 
In the tensor product the 7-matrices act on Dirac indices, while the t- 
matrices 



= (34) 
act on flavor indices. The projection operators Pjf^'^ are given by 

^f^ = ^[l®l±7^75®t5tM]- (35) 
The Dirac spinors q are obtained by the transformation 

q = q^o (36) 

whose jacobian is the function J' introduced in Eq. @. The charge conju- 
gation transformation 

q = C-'q' 

q = -C^q' (37) 

where 

C = 7o72 ® tot2, (38) 
exchanges x and y according to 

a; = 72 ® tot2 y' 

y = -72 ® tot2 x', (39) 

justifying our definition (J^) of the baryon charge. Now we reinstate Uq. The 
partition function takes the form 

Z = J[dqdq]exp{-lQSg) , (40) 

with the action 

Sq = Kj2{^w{-l^)Q{nohoPt^Uo{no)q{no + l) 

no 

-exp(/i)g(no + lhoPi^^U+{no)q{no) + q{no)N{no)q{no)} . (41) 

The factor 16 in front of the action accounts for the fact that the volume 
element with Kogut-Susskind fermions is 16 times larger than in the Wil- 
son case. The chemical potential appears in the same way as with Wilson 
fermions apart from the sign. This is due to the fact that with the con- 
ventions adopted the quarks propagate backwards in time with the Wilson 
definition and forwards with the Kogut-Susskind one. 
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5 The spin-diagonal basis 



We want to explore the possibility of using the results of the previous Section 
to construct QCD at finite density in the spin-diagonal basis. To do this the 
gauge fields must be defined on the hnks so that the block link operators 
Wf^{no) have matrix elements 

W^{no,n) = C/^(2no,2n)C/^(2no + e^,2n + e^). (42) 

The relation of the coordinates of the blocks to the coordinates of the 
sites (not to be confused with the Grassmann variables of the previous 
Sections) is 

= 2n^ + r]^, r]^ = 0, 1. (43) 
Then the action in the flavor basis becomes 

^1 = KY.{eM-I^MnohoPt^Wo{no)q{no + l) 

no 

-exp(/x)g(no + lhoPt^W+{no)qino) + q{no)N{no)qino)} (44) 



with 



-pWrj-%+(no))] + ^101 + 75® ^5^0- (45) 

But the gauge fixing Wo ~ 1 has a nontrivial Fadeev- Popov determinant and 
wc cannot say that the transfer matrix is positive definite before its effect 
is taken into proper account. Wc then proceed in the following way. We 
consider first the free case and then the interacting case to show the nature 
of the practical complication if it can be shown that transfer matrix remains 
positive definite in the gauge Wo ~ 1. 

In this Section wc do not find convenient our summation convention and 
therefore we will abandon it. 

5.1 The free case 

The spin-diagonal and the fiavour bases are related by 



qn;aAnor = ^ ErS^(2n + vY, V, = 0, 1, (46) 
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where 

rS = (7r72^^73^^7o^°).,/. (47) 

Performing this transformation in Eq. (0) with Wf^ = 1 we find the action 
in the spin-diagonal basis 



+ 5S^ I , 
(48) 



where 

a^{x) = (_i)-i+-2+...xm-i (49) 

and 5S^ is the contribution of the chemical potential. It is obtained by the 
tranformation of the corresponding contribution in the flavor basis 

5S, = 16ir^{-(l-exp(-/i))gK-l)7oPi"^g(rio) 

no 

+ (1 -exp(/i))gK)7oPo^^^g(^o - 1)} (50) 

and is 

6S^ = K {-(l-e-^)ao(x)^(a;-eo)^(x) 

+ {l - e'')ao{x)i^{x)ip{x - eo)} . (51) 

Note that the sum extends only over the sites with even temporal coordinate. 
Finally the total action in the spin-diagonal basis is 

= ao{x) [ip{x)ilj{x + eo) - ip{x + eo)tp{x) 

XQ euerijX 

+e~^tlj{x — eo)tlj{x) — e^il){x)il){x — eo)| 

+ ^ < -ft' ^ aj{x) ii{x)ip{x + Cj) - tp{x + ej)i'{xo) 

+m4{xo)ij{xo)} . (52) 

The chemical potential appears only on half the sites. We get a uniform 
coupling by the following change of variables in the sites with even temporal 
coordinate 

ip{n) ip{n) exp(/i/2), ip{n) ip{n) exp(— yu/2), no even. (53) 
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The quark action assumes then the standard form 

= ^ \^Kao{x) exp{-fi/2)ilj{x)ilj{x + eo) - exp{fi/2)ilj{x + eo)ip{x 

X 

+K'^aj{x) ip{x)%jj{x + Cj) — ip{x + ej)ip{xo) + rmlj{x)ip{x)^ . (54) 

The chemical potential now appears in all the sites in the Hasenfratz-Karsh 
form, but its value is halved. 



5.2 The interacting case 



When the gauge fields are present we must change the transformation (^61) 
to make it gauge covariant. This can be done by introducing a string of link 
variables Vri{n) connecting the site 2n to the site 2n + t] 



(55) 



A discussion of different ways to construct Vjj{n) can be found in |TT|. For 
instance we can assume 

V,in) = [U,i2n)r [f/2(2n + Si,,,e,)]'^' [Usi2n + S,,,,e^ + S^,,,e2)r 

[Uo{2n + 5i,^,ei + 6i^r^,e2 + Si^r^e^)]'" . (56) 

The contribution of the chemical potential to the action in the presence 
of the gauge fields can be read off from Eq. 



6S, = KY.{-{l-expi-fi))q{no-lhoWo{no-l)Pt'^q{no) 

no 

+ (1 - exM)q{nohoPy^W^ino - l)g(no - 1)} . (57) 
After the transformation ( ^5] ) this becomes 
6S^ = Kj2oio{r]){i'^-(^wil^))K,o 



n.ri 



■i){2n + 7])V+{n)W^{n - eo)K,+eo(^ - eo)^(2n + r/ - eo) 
-(1 -exp(-/i)) 5^0,1 

■4j{2n + r^)V+{n)Wo{n)V^.,,{n + eo)^(2n + + to)} . (58) 
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Therefore the total action in the spin-diagonal basis is 

n.ri { fi 

—B^{n, ri)ijj{2n + r] — e^)] + mijj{2n + ■r])ip{2n + ?7)| (59) 

where 

Ao{n,r]) = Sr,o,iexp{-iJ,)Wo{n)Vrj-eo{n + eo) + 5r,o,oVrj+eo{n), 
Bo{n,r]) = 5^o_oexp(/i)Vro+(n - eo)V^+eoin - eo) + 5r,o,iK?-eo('^), 
Aj{n,r]) = Sr,^^iWj{n)Vr^^e,{n + ej) +6r,^^oVr,+e,{n), 
Bj{n, 7]) = 5^.fiWj^{n - ej)Vn+e,{n - ej) + 5,,^-,iK,_e, ('^)- (60) 

Again the chemical potential appears only on half the sites and we get a 
uniform coupling by the change of variables ( |53D which changes Aq, -Bo 
A'q, B'q according to 

Ao{n,riy = e~^/2 [(5^o,iVFo(^) Veo(^ + eo) + 5„o,oK,+eo(^^)] 

Boin,!])' = e^l'^ 5^o,oW^O^(^ - eo)K7+eo('^ - eo) + 5r,o,l^'?-eo(^) • (61) 

The chemical potential now appears in all the sites in the Hasenfratz-Karsh 
form, but its value is halved. 

6 Summary 

We have seen that we can couple the chemical potential linearly to the Fermi 
fields, conserving its meaning of a Lagrange multiplier in the transfer matrix 
formalism, which is equivalent to the hamiltonian formalism in statistical 
mechanics. Starting from this definition we have derived the Hasenfratz- 
Karsh action showing that the exponentiation of the chemical potential is 
a property of the Grassmannian kernel of fermionic operators, and is not 
related to the exponential dependence of the gauge fields. 

In the case of Kogut-Susskind, the value of the chemical potential depends 
on the basis, and this can be important when the chemical potential has a 
physical interpretation, which in the present scheme is only possible in the 
flavour basis. 

The formulation of QCD at finite baryon density in the spin-diagonal 
basis cannot be considered settled. Within our approach one should prove 
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the reflection positivity of the transfer matrix in the gauge Wq ~ 1. Then 
it might also be possible to justify the minimal coupling in the spin-diagonal 
basis, by appropriately changing the gauge coupling in the flavour basis. 
Alternatively one should use a transfer matrix constructed directly in the 
spin-diagonal basis. 

We postpone to separate paper the discussion of the uniqueness of the 
path integral for a given transfer matrix. 
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